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We study the role of semi-classical QCD vacuum solutions in high energy scattering by 
considering the instanton contribution to hadronic cross sections. We propose a new type 
of instanton-induced interactions ("instanton ladder") that leads to the rising with energy 
cross section cr ~ of Regge type (the Pomeron). We argue that this interaction may 
be responsible for the structure of the soft Pomeron. The intercept A > is calculated. 
It has a non-analytic dependence on the strong coupling constant, allowing a non-singular 
continuation into the non-perturbative region. We derive the Pomeron trajectory, which 
appears to be approximately linear in some range of (negative) momentum transfer t, but 
exhibits a curvature at small t. Possible role of instantons in multiparticle production is also 
discussed. 



I. INTRODUCTION 



QCD vacuum is known to possess a rich structure. A striking example is provided by the existence of 
(Euclidean) classical solutions - instantons ||], which are responsible for non-trivial topological properties 
of the theory . While the influence of instantons on the properties of hadrons and their interactions at 
low energies have been extensively studied [^jj^,^ (for a review, see ) , the role of topological effects in high 
energy collisions is still an open and fascinating problem j8| . 

In electroweak theory, significant interest was excited by the possibility of baryon number non-conservation 
caused by instantons at the collision energy above 10 TeV ||9|-[l4|. In QCD, the instanton contribution to 
the structure functions of deep-inelastic scattering was studied in Refs . |l5|,[l^, while the cross section for 
gluon-gluon scattering mediated by instantons was considered in Refs. |13| , |17| -|19l- 

Recently, it has been proposed |2^ that the existence of semi-classical solutions in QCD can have important 
consequences also for hadron scattering at high energies, inducing a non-perturbative, and possibly dominant, 
contribution to the scattering amplitude. This approach was shown to lead to the "soft" pomeron with the 
intercept ap — 1 ~ 0.1 which has a non-trivial dependence on the strong coupling and the numbers of colors 
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and flavors 0. The method used in that work was based on low energy theorems of broken scale invariance 
p2| , and did not assume any specific form of the semi-classical solutions. This approach was based on the 
low energy theorems for the trace of the QCD energy-momentum tensor taken in the chiral limit of 
massless quarks 

^ ~:^F'^"^K0, (1) 

where b = {llNc — 2Nf)/3 and F"""^ is the gluonic field strength tensor. The low energy theorems derived 
in 1221 state that 



* J d^x{o\Te';{x)e:mo) = -4(o|^?;:(o)|o). (2) 

Comparing Eqs. ^) and (||) one can conclude that the low energy theorem (||) can only be satisfied by a 
strong gluonic fieldn 

A, - -. (3) 
9 

The strong field of Eq. ^ is usually associated with the quasi-classical solutions of QCD. This quasi-classical 
field has been employed in |2C| ] to provide gluonic interactions, which, after being iterated in the t-channel 
gave rise to a pomeron-like behavior of the corresponding cross sections. 

Here we would like to develop the idea proposed in [g^l a-nd explore the effects of quasi-classical fields on 
high energy scattering. Therefore, for the purpose of this paper we will consider a particular type of this 
classical field: we will assume that the field is given by the instanton solution . We will employ instanton 
fields, which would allow us to quantify the assumption of the quasi-classical fields being responsible for the 
pomeron-like behavior of the cross sections. Nevertheless we can not rule out the possibility of some other 
quasi-classical gluonic fields giving a significant contribution to the processes considered below. 

The instanton calculations are better justified in the electroweak theory, where the coupling is small and 
one can hope that quantum corrections to the classical instanton solution are small. Therefore the calculation 
that will be presented below can also be repeated for the electroweak interactions, where it would possibly 
also give a cross section that rises with energy 0. 

Shuryak has argued that the explicit use of instantons could provide a useful way to interpret and 
extend the results of po[] . Very recently, Shuryak and Zahed [^5| analyzed Euclidean scattering induced 
by the instantons, and analytically continued their results to Minkowski space using the correspondence 
between Euclidean angle and Minkowski rapidity ||2^. Their result is a constant cross section, which does 
not rise with energy. 

The purpose of the present paper is to investigate the effect of the instantons on the dynamics of gluon 
ladders in high energy scattering. Our approach here is complementary to the one taken in Ref. pO[ |; even 
though it is hindered by significant numerical uncertainties, it provides a deeper insight into the mechanism of 
non-perturbative effects in high energy scattering. As will become clear later, our treatment of the problem 
is different from Ref. [psf; while the authors of that paper consider a purely classical contribution to the 
scattering amplitude, we attempt at evaluating the leading quantum terms arising in perturbation theory 
built in the background of classical instanton fields. 

The proposed mechanism of the pomeron is illustrated in Fig. |^. The pomeron is being constructed 
by resummation of the ladder diagrams, similarly to the well-known BFKL pomeron |^,^^. The emission 
vertices of our pomeron ladder are given by the multiple gluon interaction vertices generated by the instantons 



^ Analogous non-perturbative behavior was also established in the scattering of color dipoles at low energy 
■^Unfortunately the applications of low energy theorems of Eq. depends very much on the way one subtracts 

the perturbative contribution to the correlation functions jiij. Lattice simulations have also encountered similar 

problems. We thank Al Mueller for clarifying this point to us. 
^We thank Larry McLerran for bringing our attention to this problem. 
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The vertices connect to two gluons in the t-channel and produce several gluons in the t-channel. We 
sum over all possible numbers of the produced t-channel gluons in each vertex loUlQ] . The t-channel lines in 
the ladder of Fig. will be taken as the usual gluon lines for most of our paper. In principle one should 
include the virtual corrections in the t-channel, which may lead to reggeization of the gluons |2^j2^. We 
will present an estimate of these virtual corrections. However the exact calculation and the answer to the 
question of whether gluon reggeization happens in our case will be addressed elsewhere [ p9[ . For most of this 
paper we will work in pure gluodynamics with no quarks. However, at the end of the paper we will present 
an estimate of the effect of inclusion of light quarks. 




FIG. 1. The structure of soft pomeron as conjectured in the paper. 

The multiple gluon vertices generated by instantons in the pomeron of Fig. |l| will be calculated in Sect. 

ni at the classical level. The vertices include a suppression factor of exp ^— , which is due to the classical 

instanton action. If the coupling constant is sufficiently small (a^ <C 1), then this suppression factor is also 
small 



exp 



27r 



< 1 



(4) 



and can be used as a parameter justifying the perturbative expansion of Fig. |l|. The scale of the running 
coupling in this case is set by the typical size of the instantons, so that as = as(l/po)- We assume that 
as(l/po) ^ 1: which is only marginally satisfied in QCD. Our ladder is resumming the leading logarithms 
of center of mass energy s, which in this case implies resummation of all powers of the parameter 



In s 



1. 



(5) 



Note that the dependence of our small parameter on is non-analytical, which will introduce the non- 
analytical dependence on Ug in the calculated value of the pomeron intercept. 

As energy in each rung of the ladder increases, quantum corrections to the tree level calculation of the 
vertices become important. The quantum corrections are known to modify the energy dependence of the 
2 n transition vertices P, p^ , p^ 14 1. The energy dependence of the 2 ^ n total (summed over n) cross 
section is given by the following formula 



(72 (s) ~ exp 



sph 



(6) 
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Here s is the center of mass energy of the system, Egph is the sphaleron energy. F is known only for the 
small values of its argument |14|. The lowest order term in F is given by the tree level calculations, and the 
quantum corrections provide higher order terms in the expansion of F [p| jic| , [l^ , [l^ . It is also known that 
for a large number of gluon legs (large n) the value of the saddle point is shifted. Thus for large number of 
produced gluons the vertices may not be exactly the same instanton-induced vertices as for small n. This 
only substantiates our argument that the instanton field may not be an exact mechanism responsible for the 
multiple gluon interactions in the ladder of Fig. |l|. There are many problems related to the determination of 
the function F{^/s/ Esph) ^^'^ ^"^^ '^^^ going to review all of them here. We will simply point out that 
the growth with energy of the cross section of Eq. (^) should eventually stop due to the unitarity constraint 



const 
a < exp 



27r 

-const — 
a.. 



(7) 



Moreover, it is possible that the discussed 2 ^ n cross section will fall off at the energies of the two incoming 
gluons much higher than the sphaleron energy, since the instanton effects may not be important in that 
kinematic region. Assuming that this is true throughout the paper we will approximate the multi gluon 
interactions by the tree level calculation for the energies below the sphaleron energy and we will just put it 
to be zero for ^/s > Esph- This approximation is rather crude and will be improved in the subsequent work 
p9{ , where several different behaviors of the 2 ^ n scattering cross section will be considered. In particular 
we will consider the case when the unitarity limit of Eq. (|^) has been reached. 




FIG. 2. Space-time picture of soft pomeron. 



Let us discuss the (Minkowskian) space-time picture of high energy scattering amplitude that we are 
going to evaluate. Since at high energies the scattering occurs over large longitudinal distances, in order 
to construct the amplitude one may need to have several instantons involved in the process (see Fig. |l|). 
At first glance, this seems to contradict the assumption of the dilute instanton gas |31| that we are going 
to use^. The resolution of this apparent problem is in the very fact that a pomeron exchange is not an 
instantaneous interaction. Pomeron exchange is best represented in the wave function picture of high energy 
interactions |pQ]. In that formalism in order to exchange a pomeron the colliding hadrons have to develop 



*A computation in a self-consistent instanton vacuum model would be a very interesting extension of this work. 
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large multi-gluon fluctuations, which, in turn, interact with each other, which corresponds to the pomeron 
exchange in the usual t-channel language. 

The space-time picture of our soft pomeron constructed in the spirit of the wave function formalism is 
shown in Fig. ^ Single pomeron exchange corresponds to a single gluon, which, after being emitted off 
the original hadron, propagates through space, inducing a chain of instanton transitions between different 
topological vacua and producing more gluons in each transition. The coherence length of the small- a; gluon 
in Fig. ^ is large, t ~ fc+//c^, where k is the gluon's momentum. Thus a high energy gluon carrying 
a large "plus" component of momentum can propagate over large longitudinal distances inducing several 
instanton transitions. Even in the dilute instanton gas approximation [pl| the gluon should therefore be able 
to generate many instanton transitions while traveling over large distances in the longitudinal direction at 
the given value of the impact parameter. 

If the process is viewed in the rest frame of one of the hadrons, then after several instanton transitions, 
the propagating gluon would interact with the target hadron at rest. In the center of mass frame the 
propagating gluon would find another gluon coming from the second hadron and would interact with it 
through an instanton transition, producing more gluons. 

Multiple pomeron exchanges in this picture would correspond to the case when the gluons produced in 
the interactions of Fig. ^in turn start propagating through the instanton gas and interacting ]3^ . However, 
a study of that interesting scenario is beyond the scope of our paper. Here we will just argue that these 
multiple pomeron exchanges are suppressed compared to the single pomeron exchange contribution, for the 
reasons described below. Of course they will become important as the center of mass energy of the collision 
becomes very high. 

The paper is organized as follows: after discussing the general issues related to instantons in Sect. II we 
will proceed to calculating the intercept of the pomeron of Fig. |l| in Sect. III. We will then calculate the 
slope and trajectory of the pomeron in Sect. IV and conclude by summarizing our results in Sect. V. 



II. GENERAL ISSUES 



In Minkowski space, instantons correspond to tunneling transitions between different topological vacua 
[^-Q . The instanton calculations therefore make sense only at energies much smaller than the energy of the 
potential barrier separating the vacua ("sphaleron" energy Esph)- 

The BPST instanton solution in the regular gauge is given by jl] 



2 l^a^v^w 



= 4=^, (8) 



and the corresponding field strength is 



The energy along the instanton path can be defined as 



E{t)^1 ld'x{G1,)\ (10) 



and its maximum at t = corresponds to the sphaleron energy (see [p3[): 



i?(r = 0) = ^l = -^. (11) 
p g^ ^asP 



To estimate this quantity, we take p ~ (700 MeV) ^, and as{l/p) — g^j^-T^ — 0.7; we get Esph — 2.4-;- 
2.6 GeV. 

At first glance, this value of Eg^h suggests that the perturbative treatment about the instanton solution 
is not applicable for hadron scattering at high energies \fs >> Esph- However, this conclusion is premature. 
Indeed, the space-time picture of high energy scattering discussed in the Introduction implies that the 
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scattering is described by the ladder-type diagrams, where the number of rungs is proportional to rapidity 
y ~ Ins, and the energy in each rung Emng is independent of the total collision energy (E'rung is typically 
on the order of a few GeV [|^). The crucial condition is therefore Emng << Egph, and it is likely to be 
satisfied (later, we will see that this is indeed the case). Therefore, the applicability of the instanton-based 
approach to high energy scattering is, at least a priori, plausible. 

Later we will find it convenient to use the singular gauge for the instanton field (the gauge in which the 
singularity of the potential is shifted to the origin); the corresponding expression is 



2 p'^Va, 



In momentum space, the instanton field is given by 



p^y 



(12) 
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(fc2)2 



1 



{kpfK^ikp) 



(13) 



III. SOFT POMERON: THE INTERCEPT 

In this section we will calculate the intercept of the soft pomeron depicted in Fig. |l| to hadronic cross 
sections and will calculate the intercept of the soft pomeron given by that equation. We assume that the 
scale of the running coupling constant is being set by the typical size of the QCD instanton po, and that the 
coupling constant at that scale is small as (1/po) ^ 1- This assumption is marginally well satisfied in QCD, 
but is crucial for the use of QCD instantons. In the recent year there have been proposed several scenarios 
of how the strong coupling constant behaves at large distances [ p5|p^ ] suggesting that it never becomes 
larger than 1. If the ideas stated in |35|j3^ ] are correct then our small coupling assumption is likely to be 
justified. Otherwise we note that the small coupling assumption is of course better justified for electroweak 
interactions psf . 

Assuming that <C 1 we will neglect all the usual perturbative QCD vertices as bringing extra powers 

of as- We will be resumming only the instanton-induced vertices, which give powers of exp ^— ^ 1 

and our model of the pomeron could be viewed as resumming powers of this parameter. 

The structure of this section is the following. We will first derive the intercept of the pomeron of Fig. |l|. 
Then we will estimate the magnitude of virtual corrections to the real emission diagrams of Fig. |l|. We wiU 
then calculate the numerical value of the pomeron intercept predicted by our model. 



A. Calculation of the intercept 

To be able to calculate the diagrams in Fig. |l| we have to construct the basic vertex employed there — 
the vertex coupling several gluons to each other through an instanton. This will be done in a way similar to 
the calculations of the baryon number violating amplitudes involving instantons in electroweak theory done 
in 1^] . To calculate the multiple gluon vertex we have to construct the Green function 

G:';;:.:X('?i'92;fci,...,fc„) = {2^)H\q, + q2-k^-...-K) (^S(gi)4(g2)^;ii(^i) ■■• ^m:(M>, 

(14) 

where ^Ji(A;) are instanton fields given by Eq. (^3|) and the brackets (. . .) imply averaging over the instanton 
sizes in the single instanton sector. The Green function of Eq. (|l^) is depicted in Fig. ^. It combines two 
space-like t-channel gluons carrying momenta qi and 52 with n final state gluons with momenta fci ... kn, 
which are taken to be on the mass shell. Thus, in Minkowski space, q^ = —9^, where q^ is the transverse 
component of the gluon's momentum, and kf — 0. We need the Green function to be in the momentum 
space representation because we are going to use it in the momentum space diagram calculations below. 
The difference between our Green function and the one usually considered in the calculations of 2 — > n 
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process in QCD and electroweak theory |g,|l5|,|9| is that the two on-mass shell gluons in the initial state of 
the electroweak process are substituted here by two t-channel gluons which are far off mass shell. 

Of course the Green function of Eq. describes the multiple gluon interactions through an instanton 
only at not very high energies. As energy increases the quantum corrections start becoming important jl^ , 
as was discussed in the Introduction. At very high energies of the incoming pair of gluons the quantum 
corrections modify the behavior of the Green function and the corresponding 2 — s- n cross section. In our 
approach we assume that at high center of mass energies of the hadron-hadron (or quarkonium-quarkonium) 
scattering the dominant contribution comes from the pomeron diagrams in Fig. |l] with a large number of 
rungs in the ladder. Thus the energy per rung, or, equivalently, per vertex in the ladder remains not too 
large, justifying the use of the Green function of Eq. (p^. We will return to quantify this issue below. 




FIG. 3. The multiple gluon vertex essential for construction of soft pomeron as described in the text. 
For a gluon on the mass shell the instanton field of Eq. (^3|) becomes 
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as 



0. 



(15) 



Substituting the field of Eq. ( |15| ) for each ki line into Eq. (14) and employing the full instanton field of 
Eq. ( |l3| ) for the q lines we obtain 



G:XV.:X (91, 92; fci, . . . , fc„) = {2Trf 5\q, +q2-h 



k„) I dpn{p)p'^ 



n 



aa-y 



ill 



2\2 



l--{qipfK2{qip) 



VbfIS (125 



('?2P)^^2('?2P) 



(16) 



where n{p) is the size distribution of instantons. We want to obtain an effective vertex for multiple pomeron 
interactions. For that we have to truncate the external legs of the Green function of Eq. (^. To amputate 
the propagators of the external gluon lines one has to also remove the numerators of these propagators. 
For that we need to know the gauge condition of the fields of Eq. ( p^ ) and Eq. (|l^). Usually the field of 
Eq. ( [T^ ) is derived using the Schwinger gauge condition a;^ = 0. However, we note that it also satisfies 
the covariant gauge condition = 0. Therefore we will be using it as a covariant (Feynman) gauge field 

and will employ it in the covariant gauge calculation of diagrams. Therefore, to truncate the Green function 
of Eq. (|l^) we have to just multiply it by q| fc^ ... k'^- The resulting amputated Green function is 



ra/T/li'.'.X^^i' ^2; fci, • ■ ■ , kn) = {2'K f5^{qi + 92 - fci 



-(47r)4 
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n 



{q2pf K2{q2p) 



(17) 



Now we are in the position to calculate the intercept of the soft pomeron in Fig. The calculation will 
be done in the Minkowski space. To perform it we have to first multiply the amputated Green function of 
Eq. ( p7| ) by the numerators of the propagators of the t-channel gluons in the light cone gauge (Weizsacker- 
Williams approximation) 

qi+ 92-' 



where the gluon qi carries a large + component of the light cone momentum and the gluon q2 carries a 
large — component. This is a standard procedure, which is described in ||37|| in greater detail. Then we 
should multiply the gluon line factors ki's in Eq. ( [l7| ) by the polarization vectors e^^{ki) and square the 
amputated Green function of Eq. (|lj) excluding the S function. We should keep in mind that the color space 
orientation of the instanton in the complex conjugate arnplitude is, in principle, different from the color 
space orientation of the instanton in the amplitude |]lT| , |l2| , p^[|l^ , [l9[| . Therefore we will have to average over 
all possible orientations, for which we will be using the approximation introduced in [T9| ]. We will integrate 
the obtained expression over the phase space of the produced gluons and sum over colors and polarizations 

A,; 



n 



Li 2(27r)= 



E 



where uji = k\ and sum over n runs from 1 to infinity. We also have to multiply everything by one of each 
of the denominators of the propagators of the t-channel gluons [0 

1 1 

and average over the colors of gluons qi and 52, which would give a factor of 1/{N^ - 1)^. We also have to 
include the symmetry factor of (n!)^. We end up with the following expression for the kernel of the integral 
equation describing the pomeron of Fig. |l| 



K{qi,q2) In- 

X 



1 



- 1)^ 



E 



d'^q6\q-Y,k^)n\ 

i=l 



d^k, 



i—1 ^ ' 



47r2 

a 



(4vr) 



9? ql 



Vaa- lla Vbp+l2l3 



9? ql 



U^'''ma'c.-<llaU'"'mb'0+<l2p 



9? 9o 



dpn{p) p^ 



l--{qipfK2{qip) 



^-l^{<lipfK2{qip) 



(18) 



with q = qi + 92- U°''^{£) is the matrix of global SU{2) rotations in the adjoint representation and it is 



responsible for rotating the instanton in the complex conjugate amplitude in the color space |11 19[|. We can 



perform the summation over polarizations and averaging over the color space orientations of the instantons 
using the approximation outlined in |19[ 



8 




/ n E [^a.... fc. e^(ki)\ k., 4:iki)\* « l^-^j « (19) 

for large n. Eq. (|l9|) can be understood in the rest frame of the system of two t-channel gluons. There 
the dominant contribution to Eq. 
state gluons, which leads to the {E'^ 
was justified in |l|,|l|,0. In writing 
factors of n. However, as could be seen from the calculations below, we do not have the precision to keep 
all the powers of n at this point. Extra powers of n may arise from different ways of averaging over the 
instanton sizes (see Eq. ( |29| ) and Eq. (|35|)). The integration in Eq. ( p^ ) is sharply peaked at the upper 
cutoff, as will be seen from the calculations presented below, which we will put to be E^pii. This allowed 
us to just put i?3p/j in Eq. (p^. However, the uncertainty in our knowledge of E'^^f^ may also introduce 
significant numerical changes, which we can not control at the moment. We are going to clarify those in 
the subsequent work |2^. We have also put U°'''{S,) = 1 for the t-channel gluon lines. As the left hand side 
of Eq. (|9|) should of course be a Lorentz invariant expression in Minkowski space it can only depend on 
boost invariant quantities. Thus we can see that for n = 1 the answer should be proportional to A:^ = and 
should, therefore, be zero. For higher values of n the left hand side of Eq. ( p^ ) may depend on the products 
of momenta of different gluons, which results in a non-zero answer. 



In evaluating the expression in Eq. (18) we have to note that the perturbative size distribution of the QCD 
instantons n(p) is divergent at large sizes p. However, the lattice data |l^] suggests that the distribution 
actually starts to fall off very steeply with increasing p after some critical value of p = pQ. This behavior 
could be expected by the following argument based on the low energy theorem of Eq. (^), which has become 
known as "6/4" problem The issue is the following: it is impossible to satisfy the low energy 

theorem of Eq. (y) with a purely classical field. The energy-momentum tensor of a classical field given by 
Eq. (|l|) has only one power of 6 in it, since the classical field has no information about renormalization and 
running coupling constant in it. Thus if one tries to satisfy Eq. (||) with this energy momentum tensor 
one would get an extra power of b on the left hand side which would not cancel. The commonly accepted 
resolution of this problem is to assume that the instantons in the instanton gas interact with each other, 
which leads to a different infrared cutoff on the integrals over p on the left and right hand sides of Eq. (^) 
p8[ . This leads us to conclude that QCD regularizes itself at large distances by cutting off the growth of 
the instanton size distribution through some non-perturbative, but not classical mechanism. This conclusion 
seems to be supported by lattice data ]l^ . For the purpose of this paper, in order to cure the problem we 
have to introduce an upper cutoff po in the p integrals in Eq. (p^. 

In our approximation the sum over polarizations and colors (Eq. (p9|)) became independent of fc^'s. There- 



fore we can perform the integration over the phase space of n gluons. The integration yields 

7^1 Li (n- l)!(n-2)! 

After some simple algebra Eq. ( p^ ) becomes 

K{qi , 92 ) In - = jirp^ tt^ } , d q dq n\ ( i 



X {N'^-'^Y t^^J -2g+ ■ (27r)3" (n-l)!(n-2)! 3" e^" (n!)2 \^ g 



9^ iti^ 12± 



Po 



dpn{p) p^ 



l--{qipfK2{qip) 



^-l^iqipfK^iqip) 



(21) 



We note again that the instanton size distribution is very sharply peaked around its' maximum given by 
Po- This statement is supported by the lattice data fl^ and by various instanton models ||3l|] . The size 
distribution for p < po is very close to the perturbative results ||l^,^ , which sharply increases near po . 
For p > Po the distribution falls off steeply Q, justifying our cutoff procedure in Eq. (^l|). Thus without 
loosing the accuracy of our calculations we can put p = po in the factors containing the confluent Bessel 
functions in Eq. (^). Now we can see that the transverse momentum integration factorizes in the expression 
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for the kernel of the equation for the pomeron-mediated forward amphtude (|2|). Each segment of the ladder 
in Fig. 1^, consisting of the square of the Green function of Fig. ^, gives a factor of 

i^(q^) ^ (l-l {q^Po? K^iq^p^)] (22) 



ql V 2 

to the rung of the ladder above it and a similar factor (with a different momentum) to the rung below it. 
Therefore each rung of the ladder brings in a square of the factor in Eq. ( p^ ) integrated over the transverse 
momentum 

hpl^ J (l-l{q±PofK2{q±po)^ ■ (23) 

Therefore, to obtain the pomeron's intercept out of the integral kernel in Eq. ( pT] ) we should substitute the 
integral Iq instead of the transverse momentum integration in Eq. (|2^). Iq is easy to evaluate numerically, 
which gives 

lo « 0.014. (24) 

Integration over in Eq. ( |2l] ) gives a factor of In ^ . Actually, the summation of the diagrams in Fig. |l] can 
be performed by the following simple integral equation 

(f>{qt,Y = ln{l/x)) = K{qi) + CK{qi) dy j K{qi) <P{qi ,y) (25) 

where the kernel of Eq. (El|) is K{q^,q2) = CK{q^) ■ K{q2) and constant C absorbs sum over n and 
integral over q^ in Eq. (plF(/)(g-L,y = ln(l/x)) denotes the sum of all diagrams of Fig. ^type, i. e., the 
pomeron-induced structure function or cross section. The initial condition is given by K{q^), which is not 
crucial for us and convenient for determination of the intercept. One can check that 

0(q^, y) = K{q^) ■ eM^soft Y) (26) 

is the solution of Eq. (E3) with 



A 



soft 



C J d^q^ K^q^) . (27) 



A much more serious problem is posed by the integration over q^ and summation over n ( both absorbed in 
C in Eq. (^5|) ) in Eq. ( pT] ) , which is potentially very dangerous if one allows the upper limit of the integration 
to be as large as the center of mass energy of the system. This is what Eq. (^l|) seems to suggest if one takes 
it at face value. However, our Eq. ( pl| ) is valid only at relatively low energies. As energy increases quantum 
corrections become important. As was shown in |lj the instanton-induced cross section for 2 — > n process 
in QCD falls off with energy for energies above certain threshold energy Eq due to quantum correction 
calculated in [^3|. Similar behavior results from the unitarity constraint derived in [l9| ]. Here we assume 
that this is also the case for our expression for the intercept. Namely we believe that at higher energies per 
rung in the ladder the quantum corrections would significantly slow down and later on completely overturn 
the growth of the expression in Eq. (|2^) with energy, as was discussed in the Introduction. These effects 
could be taken into account by putting an effective upper cutoff -Efp^j in the q^ integration in Eq. (^) . After 
performing all the integrations mentioned above Eq. ( plj) yields the following expression for the pomeron's 
intercept 



(iV|-l)2 a3 6e2 ^ [(n- 1)!]3 I 6ae' 



2 



A.o/* = T^T^f^ ^ V j—-^ I — ^ I ( / dpnip)p^- ] . (28) 







Employing the same argument about the sharpness of the distribution of n{p) p6[ as was used in performing 
the transverse momentum integral in Eq. (GW we could rewrite 



10 



Pa 



where we have defined 

No -- 

Substituting Eq. (|2|) into Eq. (|2|) wc obtain 

^hjNoptY (47r)6 pI 



A 



soft 



6 



dpn{p). 



n=2 



1)!]3 



6 a 



(29) 



(30) 



(31) 



Note that Nq has the dimensions of M ^ and, therefore, Eq. (|3l| ) is, of course, dimensionless. We can rewrite 
Eq. (^) in a more compact form with the help of a generalized hypergeometric function 



A 



soft 



(7V2-1)2 a3 



0^2 



1,1;- 



sph Po 



6 a 



1 



(32) 



One could be a little more careful in evaluating Eq. (|2^). We can use the perturbative expression for the 
size distribution of instantons 



where d is given by 



d = 



nip) = d — 

p5 



0.466e-i-679iv. 
(iVe-l)!(Wc-2)! 



2No 



Po 



1.510 



-3 



for 



A^. = 3 



(33) 



(34) 



and b = (11/3) Nc for pure gluodynamics. We also used 1/po as a renormalization scale in Eq. (p3|). With 
the distribution of Eq. (p3|) we obtain 



dpn{p) p^ 



^di^-l 



2 We 



2n-4 

Po 

2n + 6 - 4 ' 



(35) 



which, after plugging it into Eq. (|22) and summation over n in it yields the following expression for the soft 
pomeron's intercept 



^soft 



nlod'^ 



(4^)6 1 E^^pI 
81 6e2 



9 9 11 n nE^phPo\ _ 
2' 2' ' ' 2 ' 2 ' 6ae2 



(36) 



where we have explicitly used & = 11 for A^c = 3. The perturbative instanton size distribution of Eq. ( p3| ) 
has been shown to fit very well the recent lattice data for p < po Therefore the estimate of the 

intercept given in Eq. (|3^ ) is probably more precise than the estimate from Eq. (p^. Though in both cases 
we have to make the same assumptions to simplify the transverse momentum integration in Eq. (|T]). We 
also had to put a cutoff on p integration in arriving to both Eqs. (|3^ ) and (^6|), which is not a very bad 
approximation judging from lattice data ||l6|], but still is an approximation. In Sect. IIIC we will compare 
numerical predictions of Eq. ( ^ ) and Eq. (|36[) for the soft pomeron's intercept. 

It should be stressed that the approach developed here in Eq. (^) - Eq. ( |36| ) is an example of the multipe- 
ripheral approach (see Ref . and references therein) . It gives all typical properties of the multiperipheral 
"ladder" diagrams, namely, the typical transverse momentum of produced particles (gluons) which does not 
depend on the total energy {q'^po « 2 -;- 3 in our case, which follows from Eq. (|2^)) and the average energy 
s between the produced two bunches of gluons in Fig. ^ which also does not depend on the total energy and 
is equal to 



^soft 



\n{s/Ef 



sph 



1 



(37) 



These features of our approach cannot be spoiled by virtual corrections but, of course, they have to be taken 
into account when one considers the numerical value of the pomeron intercept. 
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B. Virtual corrections 



There are two types of virtual contributions to the soft pomeron we are discussing. One is the additional 
gluon exchanges between instantons in t-channel, an example of which is shown in Fig. ^B. These corrections 
do not introduce extra logarithms of energy (In i) in the problem. The other type of virtual contributions 
could come from reggeization of the gluon propagator due to the instanton interactions. This effects of course 
may introduce additional powers of In i . 

Let us first concentrate on the corrections brought in by the t-channel gluons as demonstrated in Fig. 
We are going to compare the contribution of the real diagram which is included in our model of the soft 
pomeron above, given by Fig. WA to the contribution of the same diagram with a virtual correction added 
to it, which is depicted in Fig. ^B. 

The difference between the graphs in Figs. and ^ is just in extra gluon line, which introduces a loop 
integral. After some simple algebra one can see that the integral brings in a factor 

K ^ (^X ( ^ [1 - h {kpo f K^jkpo)] ' [1 - HIp - k\por mp - k\po)] ' 

\ 9 J J (27r)4 (fc2)2[(p_fc)2]2 ' ^ ) 

where the labeling of the momenta is explained in Fig. |^. One of the t-channel virtual lines carries the 
4-momentum k and the other one carries p — k. In Eq. (|38| ) we included only the factors which make 
the contribution of the diagram in Fig. ^ different from the contribution of the graph in Fig. We 
also assumed that the instantons in the amplitude and in the complex conjugate amplitude have the same 
orientation in the color space. This assumption is not going to change our final answer by a large factor, 
and, therefore is good for the estimate we are going to perform here. 




B 



FIG. 4. (A) A "real" diagram contributing to the pomeron intercept calculated in Sect. IIIA; (B) An example of 
a virtual correction. 



The value of the integral in Eq. (jSg) is easy to estimate. First we note that the expression in Eq. (pSD is 
a decreasing function of p. Therefore one can obtain an upper bound on the value of the integral (pq) by 
estimating it numerically for p = 0. The result reads 

K = 0.0052;^ pi (39) 
2 a 



We have to compare the factor given by Eq. (|39|) to the factor given by the p-line in Fig. ^A, which is 
just l/(p2)2. For the estimate we can assume that approximately p « l/po (see Eq. (p^)). Therefore we 
conclude that the diagram in Fig. |^B is different from the diagram in Fig. 0A by the factor of 0.0052 7r/2 a. 
For the values of a taken at the scale corresponding to the instanton size (q;(1/po) ~ 0.6) this number is 
still very small, of the order of 0.01. We can see that resummation of t-channel virtual corrections, which 
are not enhanced by logarithms of energy can not give us a large contribution and could be neglected. 
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FIG. 5. Virtual diagram of Fig. 4B, which contribution is calculated in the text. 



One might think that the diagram of Fig. ^ taken without the s-channel giuon lines would give us a 
contribution to the mass of a gluon, thereby creating a non-zero gluon's mass. However, in constructing 
Eq. (^8|) we did not impose the condition that the gluons forming the gluon loop should be in the color octet 
state, as one should do in calculating the correction for the gluon's mass. That condition was not used since 
with the s-channel gluons present in Fig. || the color of the pair of t-channel gluons is not fixed. Now we 
have to point out that for p = the contribution of this diagram without the s-channel gluons is zero, which 
could be seen after projecting the t-channel gluons onto the color octet state. 

The second type of virtual corrections are enhanced by the factors of In i and are similar to gluon reggeiza- 
tion corrections for the BFKL pomeron |2^j2^. To include these corrections which bring in the evolution 
in In i one has to calculate the intercept of the same pomeron ladder of the type shown in Fig. |l| with two 
t-channel gluons forming a color octet combination and with the non-zero momentum transfer in the ladder. 
The two gluons in the t-channel should in turn include all the virtual corrections in them. A careful inclusion 
of the gluon reggeization proved to be a very complicated task and will be completed elsewhere p9| ]. here 
we are going to estimate these corrections. 

Let us calculate the contribution of a single rung in the ladder of Fig. |], with the ladder taken in the 
octet state with non-zero momentum transfer t = —q^ (see also Fig. |^). A direct repetition of calculations 
in section IIIA leads to the following expression for the octet channel integral Is,{q^) which replaces integral 
Iq in Eq. ( ^2|) and Eq. (|3^ ) for the intercept of the Pomeron. 



cfk. 



l(fc-lg) X (fc+^g)P 

[(A + k)T[(fc-k)T 



P0K2 [ \k+-q\po 



1 - 



2- 



2-" 



(40) 



We will now assume that the virtual corrections to our pomeron lead to reggeization of the t-channel gluons. 
Moreover, to estimate the trajectory we will assume that the expression for the gluon Regge trajectory has 
a form 



aa{q') = 1 + A.„,, . 

^0 

Eq. ( ^ ) is our estimate of the gluon's trajectory. One can see that If,{q^) ^ at ^ 0. 
To simplify the expression for the gluon's Regge trajectory we note that approximately 



l--z^K2{z) 



16 



(41) 



(42) 



Using Eq. (E^) we can easily calculate function Isiq"^), namely, it turns out that 
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256 



Ei 



2Ei 



(43) 



Eq. ( ^ ) should be substituted into Eq. to obtain our approximation for the gluon's Regge trajectory. 

Now we have to determine whether the reggeization corrections can significantly influence the value of the 
our pomeron's intercept. It turns out that gluon reggeization can change the value of A^oft since it leads to 
an additional factor in Eq. (E3), which reduces this equation to the form 
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1 



-2[aG(9i)-l] In(s7-B' J 



(44) 



where s is the energy in one rung of the ladder (see above) . As we have discussed the dominant contribution 
in the integral of Eq. (H) comes from q^pl « 2 - 3 (see Eq. (|2|)) where 2 [aG{qi = ^1 pi) - 1] ln(s/£;2p,J 
0.4 0.5, since ln(s/ii;^p^) 1/Aso/f ( see Eq. (|3^) ). Therefore, the additional form factor in Eq. (||) 
due to gluon reggeization can suppress the value of /qi and, consequently, the value of IS.aojt by a factor of 2 
3. Explicit numerical calculations confirm these expectations. A more careful analysis of the contribution 
of virtual corrections will be done later p9| ]. In the numerical estimates which will be done below we have 
to keep in mind that virtual corrections can introduce a factor of 2 3 uncertainty in the value of the 
pomeron's intercept. 



C. Value of the intercept 



The main purpose of this paper is to provide a qualitative mechanism which could account for the effect of 
the soft pomeron. This was done in Sect. IIIA. Nevertheless, we could try to estimate the numerical value of 
the soft pomeron intercept obtained in Eqs. ( ^^ and (^6|) and check how it compares with the experimental 
value of 0.08 0. 



0.08 



0.06 



0.04 



0.02 




1.2 a, 



FIG. 6. Intercept of our soft pomeron as a function of for E^-ph ~ 2.4 GeV, po = 0.3 fm and the virtual 
correctionsgiving a factor of <5 = 0.31 suppression. The horizontal line corresponds to the phenomenological pomeron 
intercept 



At the same time one has to be very careful with the numerical estimates of the intercept in Eq. ( |32[) and 
Eq. (^6|). It is important to realize that the intercept of Eq. (^2|) strongly depends on the position of the 
maximum of the instanton size distribution po and on the value of the integral of this distribution over all 
/o's given by Nq. The value of the intercept in Eq. ( |3^ ) depends on the parameter po since it determines the 
scale of the strong coupling constant. The current knowledge of these parameters is poor [ p^|jl7| , |4^ , ^ . The 
value of Nq seems to be undetermined up to an order of magnitude [E2| . Uncertainty in the value of po brings 
in an uncertainty in the value of the strong coupling constant a{lJpQ). Finally a very important question 
concerns the value of the cutoff, which was taken to be Esph in Eqs. ( p^ ) and (^6|). The cutoff corresponds 
the maximum invariant mass of the particles produced in each vertex. In electroweak physics it correspond 
to the energies when the growth of the cross section of the 2 ^ n instanton-induced process slows down and. 
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maybe, even starts decaying. In QCD it may correspond to the energy at which the unitarity hmit of the 
2 — > n process is reached [M . The exact value of the energy at which this happens is still not determined 
precisely |l^Jll|,|l^,|l^,^ . It is one of the most crucial assumptions of this work that the turnover of the 
cross section does happen in QCD and that the corresponding invariant mass i^fp^ is not too large. 

The intercept of the pomeron calculated above (Eq. (^)) has been derived for pure gluodynamics. If we 
want to compare our result to the experimentally measured one p[ | we have to include quarks in the theory. 
They would introduce suppression of the amplitudes through the Faddeev-Popov determinant [|| . Quark 
lines could be included in the scattering in two ways: they could be produced from an instanton into the 
final state, in which case the suppression factor they introduce is small, of the order of 10~^, so this type of 
quark effects could be neglected. Alternatively the quarks could be absorbed in the gluon condensate Q]. 
In the latter case the suppression factor in the distribution of instantons Eq. (js^) is Q] 

l.'i(m,p-'^-^{Q\qq\Q)p^\. (45) 



n 



3 



The suppression for the intercept of Eq. ( ^ ) is given by the square of the factor in Eq. (^). For our estimate 
we will just put p — pqih Eq. ( ^5| ) and use (0|5g|0) = — (0.25GeV)^. Taking the product over three flavors 
we obtain a suppression factor of 0.023 for the intercept. 

Keeping in mind all of the above mentioned limitations we substitute the value of Egph — 2.4 GeV and 
Pq = O.^fm from into Eq. ( |36| ) and plot the resulting value of the pomeron's intercept as a function 
of the coupling constant a^. We also included suppression introduced by virtual corrections by including 
an extra factor S on the right hand side of Eq. (|3^). The resulting pomeron's intercept is related to the 
intercept of Eq. (E^) by the following relation 



n 1-3 [m,po-^(0|5g|0)p3 



A.o/t (46) 



In the plot shown in Fig. || we have put i5 = 0.31. Of course we have to admit that the plot in Fig. ^depends 
on the value of the suppression factor 5. It is also very sensitive to the cutoff Esph ■ A work on putting some 
constraints on these parameters is under way With the progress in the field of instanton models the 
exact numerical value of the soft pomeron intercept might be calculated with more confidence and precision. 

The maximum of the intercept depicted in Fig. ^ is close to the experimental result of 0.08, shown there 
by a horizontal line. The value of the coupling constant at the scale set by the instanton size po is near 
the point on the plot of Fig. ^ where the maximum is achieved. Thus our value of the intercept should be 
about 0.08. Also, the pomeron's intercept seems to vanish at small distances (large momenta), in agreement 
with experimental data pl[ . Also our model of soft pomeron allows a non-singular continuation of the value 
of the intercept into the region of large distances and large coupling constant. Even if the strong coupling 
constant becomes very large at small momenta (large distances) the intercept of our pomeron still remains 
small (see Fig. ^). That way the strong interactions are able to yield us with a small pomeron's intercept 
even in the kinematic region where they are very strong. Here we have to recall that strictly speaking our 
quasi-classical approach is applicable only when the coupling constant is small. Nevertheless, the smallness 
of the intercept at large coupling is an interesting result, supported by phenomenological observations p| . 



IV. SOFT POMERON: THE SLOPE AND TRAJECTORY 



In this Section we are going to calculate the slope and the intercept of the soft pomeron conjectured above. 
We will include contribution of the virtual corrections by putting a factor S in the expression for pomeron's 
trajectory (see Sect. IIIB). The quark contribution will also be treated similarly to Sect. IIIC. Therefore the 
slope and trajectory of the soft pomeron which will be derived below will result only from the real emissions 
part of its kernel. This will correspond to summing the diagrams of the type shown in Fig. |l|. 

To calculate the slope of the pomeron considered here we will follow the standard procedure. Let us 
consider a pomeron with non-zero momentum transfer described by Mandelstam variable t. A rung of the 
off-forward pomeron's ladder is shown in Fig. |^. One t-channel gluon carries a momentum k + ^q, whereas 
the other one carries momentum k — ^q. Thus the momentum transfer is t = —q^. 



15 



FIG. 7. Soft pomeron with non-zero momentum transfer. 



It is easy to see that the only difference between the pomeron ladder with non-zero momentum transfer 
and the t — pomeron considered in Sect. Ill is in the transverse momentum integral. An easy calculation 
shows that for off-forward pomeron the integral Iq in the expression for the forward pomeron intercept ( ^ ) 
should be replaced by 



1-^ (k+ll) pIk2 \\k+\q\pa 



P0K2 (\k-lq\Po 



(47) 



which naturally reduces to Eq. (g3|) when q = 0. To obtain the slope of our pomeron we need to expand 
Eq. (^) up to the quadratic term in q"^ . After expansion of Eq. (^) in q^ and integration over the angles 
we obtain 



/>oo 

I = lo + Trg^po / dkpn L{kpo) 
Jo 



(48) 



where 
L(fcpo) 



^2(kpor ~ ^(^'^o)'^2(fcpo) ) [32 + (fcpo)'i^i(fcpo)' + 4(/cpo)'i^2(fcpo) + 8(fcpo)'i^|(fcpo) 



-2ikpofKiikpo) + 2 {kpo)*Ko{kpo) ( 1 - ^(kpo fK^ikpo) ) - 20{kpofK^{kpa) + 2{kpaf K2{kpo)K^{kp^) 



-{kpofKlikpo) + 2ikpofKi{kpo) (-10 + {kpofK2{kpo) + {kpof K^ikpo)) + 2{kpo)^K^{kpo) 

~ikpofK2{kpo)Ki{kpo)] . 



(49) 



The integral in Eq. (^8|) is dominated by the small values of the argument, as could be seen from the behavior 
of the function L{kpo). We can expand the function L{kpo) of Eq. ( ^ ) to obtain 



L(fcpo) ~ - 



1 

256 k po 



(50) 
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when the argument is small. The full integral of Eq. (^7|) is convergent, whereas the integral in Eq. ( [48[ ) 
seems to be divergent in the infrared. Therefore to regularize it after we plug Eq. (^0|) in it we have to use 
g as a lower cutoff and 1/po as an upper cutoff. The result yields 



""221 1 

— g>5 In . 

256 qpQ 



Recalling that a pomeron trajectory is given in terms of its' intercept and slope by 

A(i) = A + a'i 

and comparing ([54) to Eq. (Eq) and Eq. (pit) we can conclude that for our pomeron 



a' 
A 



-Pi 1^.1 



256/o qpa 



(51) 



(52) 



(53) 



Since there is still some q dependence left in Eq. ( p3[ ) we conclude that our pomeron's trajectory is not quite 
linear in t and the notion of pomeron's slope is not very well defined for it. A similar result was obtained for 
the soft pomeron trajectory in |^5|] using Regge theory, which seemed to confirm experimental data of ||4^ . 
The conclusion of |Q was drawn out of two-pion exchange model of hadronic interactions. We can obtain a 
numerical estimate of the pomeron's slope of Eq. js^ ) by assuming that In ^ ~ 1. Substituting the values 
of Jo from Eq. (^) and po = 0.3 /m = l.bGeV^^ into Eq. (H) we obtain 



A 



2.0Gey-^ 



which is remarkably close to the experimental result for the soft pomeron 



a' 
A 



exp 



0.25 



(54) 



(55) 



8D (or either Eq. (|3g) or Eq. {^) can give 
could be estimated numerically for positive 



In principle the integral of Eq. (^ ) substituted in the Eq 
us the full trajectory of our pomeron. The integral in Eq. (^ 
values of q^, which corresponds to the negative values oi t — ~q^. However, for positive t the integral in 
Eq. (^) has a singularity and if taken literally would be simply divergent. The resolution of this problem 
is the following: one should first perform the integration in Eq. ( ^ ) for negative t and then analytically 
continue the results into the region of positive t. This should be done analytically and the exact analytical 
treatment of the integral in Eq. (^7|) seems to be rather complicated. Instead we will use an approximation, 
which preserves all the main qualitative features features of the answer. 

Using Eq. ( ^2|) in Eq. ( ^7|) and performing the integration over k yields the following expression 



lit) 



256 



2 



2Ei 



2 



Ei 



(56) 



We can check the quality of the approximation done to obtain Eq. 
with Jo in Eq. E^. Eq. ^ gives 



( 15^ ) by comparing its' value for t — Q 



7(0) « — 
^ ' 256 



0.012, 



(57) 



which is very close to the exact value Iq = 0.014 in Eq. (p^. If we expand Eq. ( |5^ ) for small t up to the 
terms linear in t we would exactly recover the expression for the slope of the trajectory given by Eq. (]53|)! 
Thus Eq. ( |5^ ) provides us with a very good approximation of the integral in Eq. (^7|). 

Our soft pomeron's trajectory could be obtained by substituting I{t) from Eq. ( |56| ) into Eq. ( |3^ ) instead 
of /q. This gives 



A.o/t(t) 



{N^ - 1)2 
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_^ (4^)6 1 Elpf^pl 
81 6e2 
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(58) 



After taking into consideration the virtual corrections and the quark contributions the answer for the 
pomeron's trajectory becomes 



Ap{t) = S 



U 1.3 i^mgPo^'^iO\qq\0)pl 

q—u,d,s,... 



Aso/t(i), 



(59) 



where Asoftit) is given by Eq. ( |58[ ) 

The soft pomeron's trajectory of Eq. ( |59| ) is depicted in Fig. |[ It is plotted for pq = 0.3/m, Egph — 
2AGeV, 5 — 0.31 and a ~ 0.75, i.e., the same values as were used for the estimates of the pomeron's 
intercept in Sect. IIIC. 
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FIG. 8. Soft pomeron's trajectory, t is measured in the units of pg- 

As could be seen from Fig. |^ the soft pomeron's trajectory never becomes negative for negative t, which 
means that the total cross section generated by the soft pomeron exchange is always growing with energy. 
For positive values of t the trajectory first starts growing but then slows down and even turns over and falls 
off. It does not reproduce the linear growth behavior which would be predicted by Regge theory. However, 
the linear behavior of the trajectory for positive t is usually associated with confinement in QCD. It is a 
well known fact that the instanton-induced effects can not account for QCD confinement |^ . Therefore one 
should not expect that the pomeron constructed out of instantons should exhibit such confinement features 
as a linear trajectory for positive t. It is most likely that different non-perturbative physical mechanisms 
play an important role in that region. 



V. SUMMARY AND DISCUSSION 



In our approach we envision the soft pomeron as a chain of topological transitions induced by high 
momentum gluons. This is the main point of the paper. The corresponding pomeron ladder is depicted in 
Fig. |l| The multi-gluon vertices of the ladder are generated through a strong quasi-classical vacuum field. 

To quantify the pomeron illustrated in Fig. |l] we have used the instanton fields, which have an explicit 
analytical form which we could employ. The resulting expressions for the pomeron intercept are given 
in Eqs. ( p^ ) and (p6|), corresponding to different ways of averaging over the instanton sizes. We have 
estimated the virtual corrections and showed that they are likely to change the value of the soft pomeron's 
intercept by a factor of 2 3 at the most. The numerical value of the intercept has been plotted in Fig. 
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^ as a function of the running coupling constant. The resulting intercept at as{l/ po) — 0.7 is very close 
to the phcnomenological value of 0.08 The intercept falls off at small distances, which agrees with 

phenomenological observations of At short distances (small momenta) the dominant contribution to 
the structure functions comes from the hard pomeron pl| . The work on understanding how our soft pomeron 
merges with hard (BFKL) pomeron in that region is under way |^7| . At large distances (small momenta) 
the pomeron's intercept also decreases. That way, in the region where the coupling constant is large and 
the interactions are very strong the pomeron's intercept can still be small. Thus our results could provide 
the resolution of the long-standing puzzle of the smallness of the soft pomeron intercept in the strong 
coupling regime. 

Finally, we have calculated the slope of the soft pomeron in Eq. (M) . The value of the slope turned out to 
be close to the phenomenologically suggested value of 0.25 GeV^^ The full trajectory of the pomeron 
has also been derived (see Eq. (^) and Eq. (|59|)) and plotted in Fig.^P Since the instantons can not account 
for confinement in QCD we do not expect a linear trajectory for positive t. For negative as the absolute 
value of t is getting larger the trajectory ceases to be linear and exhibits some curvature. 

We can also estimate the multiplicity of the produced particles in a hadron-hadron scattering event 
mediated by our pomeron of Fig. |l|. The particles would be produced by multipcripheral mechanism. From 
the estimate of Eq. ( p7[ ) we can conclude that the density of instantons per unit of rapidity is rather small. 
However, since there are several particles emitted coherently from each instanton, we may expect correlations 
in particle production. 

Let us estimate the typical number of particles produced on an instanton in the ladder of Fig. 0. For that 
purpose we will consider the series of Eq. ( ^ ) and estimate the average n in that series. That would be a 
typical multiplicity of gluons produced off an instanton. A simple calculation yields 

{n) = 1 + -r^^ ) J - 60 

oF2(-;i,i;^) 

For the values of the parameter we have used above in estimating the pomeron's intercept, Egph = 2.4 GeV, 
po = 0.3/to, as « 0.7 (see Fig. we obtain (n) « 3. This result agrees with the estimate of pOj . 

While our numerical results are very sensitive to a number of instanton-related parameters, we believe 
that our approach can provide an interesting new insight in the mechanism of high-energy scattering and 
particle production in QCD. Further theoretical and experimental developments are clearly needed to clarify 
the relationship between the soft pomeron and the properties of the QCD vacuum. 
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